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Abstract 

Motivated by recent studies of deep inelastic scattering (DIS) off the Af = 4 super- Yang-Mills (SYM) 
plasma, holographically dual to AdS§ x S 5 black hole, we in this note use the spacelike flavor current 
to probe the internal structure of one holographic quark-gluon plasma, which is described by the Sakai- 
Sugimoto model at high temperature phase (i.e., the chiral symmetric phase). The plasma structure function 
is extracted from the retarded flavor current-current correlator. Our main aim here is to explore the effect of 
non-conformality on these physical quantities. As usual, our study is under the supcrgravity approximation 
and the limit of large color number. Although the Sakai-Sugimoto model is non-conformal, which makes 
the calculations more involved than the well-studied A/"=4 SYM case, the result seems to indicate that the 
non-conformality has little essential effect on the physical picture of the internal structure of holographic 
plasma, which is consistent with the intuition from the asymptotic freedom of QCD at high energy. While 
the physical picture underlying our investigation is same as the DIS off the A/"=4 SYM plasma with(out) 
flavor, the plasma structure functions are quantitatively different, especially their scaling dependence on the 
temperature, which can be recognized as model-dependent. As a comparison, we also do the same analysis 
for the non-critical version of the Sakai-Sugimoto model which is conformal in the sense that it has a constant 
dilaton vacuum. The result for this non-critical model is much similar to the conformal J\f = 4 SYM plasma. 
We therefore attribute the above difference to the effect of non-conformality of the Sakai-Sugimoto model. 
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I. INTRODUCTION 



In heavy ion collision, which is now experimentally studied at RHIC and LHC, the so-called 
quark-gluon plasma seems to be strongly interacting and behaves like perfect liquid (![, which is 
greatly different from the previously recognized weakly-coupled one. This brings non-perturbative 
investigations of hadronic matter at high temperature and high density produced in heavy ion col- 
lision into an urgent stage. Though the lattice method can calculate some properties of strongly- 
coupled system, it is still constrained to extract some static quantities such as hadron mass spec- 
trum and thermodynamical behavior. What is worse is that when adding finite density or chemical 
potential into the thermal QCD, the lattice calculation usually confronts the notorious sign prob- 
lem. Therefore, improvement in theoretical understanding of strongly-coupled quark-gluon plasma 
(sQGP) should not only go beyond the traditional perturbative QCD (pQCD) approach but also 
reveal some properties out of equilibrium, such as transport properties, dispersion relation, high- 
energy scattering and so on. 

Gauge/gravity duality Q states that strong coupling gauge theory can be mapped to weak 
coupling gravity with a negative cosmological constant in the limit of large 't Hooft coupling and 
large N c , where N c is the number of color of gauge theory. Although the gravity dual of the 
realistic QCD has not been established till now, one expects this duality to be of great importance 
in understanding some non-perturbative properties of QCD or at least some universal features of 
strongly-coupled system. In fact, over the last decade, using gauge/gravity duality technique to 
study properties of sQGP has gained great success Q. However, most of the studies have focused 
on static or hydrodynamic properties at large-scale or long-time compared to the inversion of the 
temperature of the system. Therefore, it is of great interest to study hard probe of the plasma 
and reveal its internal structure, which should be in contrast with the parton picture of a single 
hadron in pQCD. Studies on this topicQ was originally proposed in || for 7£-current scattering off 
a dilaton hadron in hard-wall model and later generalized to the plasma case without flavor in |?J 
and with flavor in ||. 

The main lessonP] from these investigations is that at high energy the 7£-current or flavor current 
probes the partonic behavior of the plasma, giving nonvanishing plasma structure function; while 
at low energy the current is not absorbed by the plasma, indicating vanishing contribution to the 
plasma structure function^]. Besides the plasma structure function, another important quantity 

1 For studies on deep inelastic scattering under the technique of gauge/gravity duality, see, e.g., refs Q and ||. 

2 For details on discussions of the structure function and the partonic picture of holographic quark-gluon plasma, 

see §-|. 

3 When taking into account the non-perturbative tunneling effect of current when encountering a narrow and high 
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is the so-called plasma saturation momentum, which is the critical energy defining the transition 
from weakly quasi-elastic scattering to high-energy deep inelastic scattering. In other words, the 
partonic picture for holographic plasma exists only when all partons have transverse momenta 
below the saturation momentum. 

Since previous studies are mainly focusing on the D3-brane geometry, which is dual to the 
conformal jV=4 SYM, we in this note^ use one non-conformal gravity dual model of QCD to 
probe the effect of the non-conformality on the plasma structure. To be specific, the model under 
consideration here is the transversely- intersecting D4/D8/D8 brane system, which is now usually 
referred to as the Sakai-Sugimoto model [10|. It is one of the most successful holographic QCD 
models from the top-down approach in realizing some phenomena of low energy QCD such as 
confinement /deconfinement phase transition^, non-Abelian chiral symmetry breaking and vector 
meson dominance. The temperature under gauge/gravity duality approach can be realized by 
extending the color brane geometry to a black hole [12]. The high temperature phase means that 
the chiral symmetry is restoring, which is denoted by parallel profile of flavor D8/D8 branes. Note 
that the temperature will be much smaller than the four-momentum of the flavor current, since in 
this note we are focusing on high energy scattering to probe the internal structure of holographic 
plasma. The setup in this note is similar to || , but due to the fact that the induced metric on the 
flavor worldvolume reduces to AdS^ x S 4 black hole, it is essential the same as that of the 7£-current 
DIS off the JV=4 SYM plasma. Thus, the procedure here is similar to [fjj and much simpler than ||. 
To extract the structure function, we need to study the flavor current propagation in holographic 
plasma. According to gauge/gravity duality, this can be achieved by studying the U(l) flavor 
gauge field in curved spacetime. The corresponding action and equation of motion determining 
this dynamics are encoded in a Maxwell theory in curved five dimensional spacetime. On the 
other hand, the parton structure function has standard definition in quantum field theory, which 
is encoded in the retarded current-current correlator. One main task of gauge/gravity duality is to 
calculate the retarded Green function for finite temperature field theory from dual gravity following 



the recipe in [13]. In this sense, the information of the structure function is totally encoded in 
the solution of the above-mentioned Maxwell equation in five dimensional curved spacetime of 
asymptotic AdS type. One additional key point under this prescription is that the solution should 



potential barrier, an exponentially suppressed structure function can be obtained. 

4 As in the literature, we here for simplification do not consider charge density which can be modeled by time 

component of flavor gauge field in present construction and just leave this problem for further investigation. 

5 Recently, there appears one paper jllj which proves that the previous interpretation of transition between the AdS 

soliton and the black D4 brane as the strong coupling continuation of the confinement/deconfinement transition 
in four dimensional Yang-Mills theory is not valid. The authors there proposed an alternative gravity dual of the 
confinement/deconfinement transition. For details ononis topic, see the original work jyj. 



obey incoming wave boundary condition at the horizon, reflecting the full-absorption characteristic 
of black hole. Once the equation of motion for the flavor U(l) gauge field is turned into Schrodinger 
type, it will be found that their behavior is very similar to that of the A/"=4 SYM case. This is 
natural since the setup of holographic models is very general. This may also be recognized as one 
universal characteristic of gauge/gravity duality. 

Another motivation for our study in this note is to probe some universal features of the structure 
of holographic plasma, described by the intersecting D-brane system in conformal or non-conformal 
case. We have learned that the procedure and the physical picture for these calculations are 
universal, which is in part due to the unified approach of gauge/gravity duality in dealing with 
strongly-coupled problems. The specific form of the structure function is model-dependent, which 
may allow us to resort to the experiments to judge which model is much closer to physical reality. 
For the two structure functions, we find that in the non-conformal background they also satisfy 
the Callan-Gross relation in the limit of large Bjorken variable being defined later, which is same 
as in the AA=4 SYM case. Therefore, it is reasonable to conjecture that this relation should be 
universal for holographic plasma described by the intersecting D-brane system, having nothing to 
do with the conformality of holographic background. 

Since the flavor gauge field we are focusing on in this note can also be considered as the 
gravity dual of vector meson, our study can then be regarded as the completion of previous studies 
of mesonic quasinormal modes in [[Df ]. Here, we extend these studies to high momentum and 
high frequency limit, in contrast to previous hydrodynamic behavior or just high frequency limit. 
However, here we will not go into the detailed numeric extraction of mesonic quasinormal frequency 
in high frequency and high momentum limit, and leave this task for future direction. 

We will also do the same analysis for the non-critical version of the Sakai-Sugimoto model for 



comparison. This non-critical model was proposed in [15| to overcome some drawbacks of the 
critical Sakai-Sugimoto one, which are general in critical string theory. This model can be thought 
of as a conformal one, so we expect the related results in this model should be more similar to the 
AA=4 SYM case. This expectation has been confirmed in the fact that the structure function for 
non-critical model has the same scaling dependence on temperature as that of jV=4 SYM case. 
This should be another difference between the Sakai-Sugimoto model and its non-critical version. 

Another striking point is that the stringy imprints have appeared in the final results of the struc- 
ture functions for both models considered in this note. Seemingly, this would sway the successful 
aspects of the Sakai-Sugimoto model. However, if we recall that the holographic plasma is quite 
different from the realistic sQGP, we should just take these unsatisfactory features as non-universal 



ones and focus on some universal features emergent from different holographic models. 

The rest of this note is organized as follows. In Sec. [l| we first give a brief overview on the 
Sakai-Sugimoto model as well as its non-critical version, and then give the basic equations for later 



calculations. In Sec. Ill we state detailed extraction of the structure function for the Sakai-Sugimoto 



model and list the main results for its non-critical version. Then we have a brief discussion about 



the results. Sec. IV is devoted to a short summary and some open questions. 



II. OVERVIEW OF MODELS AND BASIC EQUATIONS 

In this section we first recapitulate the Sakai-Sugimoto model and its non-critical version p0|, 



15, 16]. Then we turn to the definitions for the plasma structure function in terms of physical 
quantity in standard field theory. We will also state basic equations determining the flavor current 
propagation in these plasmas from the viewpoint of dual gravity, which are essential for later 
extraction of the structure function. We will follow the notation conventions in |jl H]. 

A. Sakai-Sugimoto model versus its non-critical version 

The bulk background geometry of the Sakai-Sugimoto model is given by a ten-dimensional 
supergravity description of N c coincident D4-branes in type-IIA superstring theory compactified 
on a circle. According to ]l2|], there are two different metrics for this supergravity, representing 
two different phases of holographic QCD. The transition between these two phases is interpreted 
as the deconfinement phase transition. Here, we just focus on the high temperature deconfining 
phase, described by following backgrounds^ 

3 

ds 2 = (|) f (-f(u)dt 2 + dx 2 + dxj) + (?\ 5 (JjfL + u * d ni\ , (1) 

3 ^ 

e ' = Mi*) ' f ^ = 1 -{-) ' t^t + P = t + - T . (2) 

The temperature of holographic plasma dual to the above background is 

i 

T=\ = ^- (3) 
P 4vri?2 

The curvature radius R of the background is related to the string coupling g s and string length l s 



via 



R 6 = 7rg s N c li = TcXlf. (4) 



6 Due to the periodic identification of time coordinate, it should be the Euclidean time. However, we here use the 
Minkowskian signature for later convenience of extraction of Minkowski space retarded Green function. 
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Here in the second equality we have defined the 't Hooft coupling constant A from the dual gravity 
side as X = g s N c . 

The above gravity background is dual to the gluon sector, and the quark sector can be introduced 
in quenched approximation by adding Nf pairs of D8 and D8 flavor branes to the above geometry 
and make them transverse to the circle along X4. In the quenched limit, Nf <C N c , the backreaction 
of the flavor branes on the background geometry can be neglected. Dynamics of the flavor sector is 
encoded in the Dirac-Born-Infeld (DBI) plus the Chern-Simons (CS) actions for the flavor branes in 
above background. However, the CS-term will be exactly zero in this note as there is no background 
for the U(l) gauge field on the flavor branes. 

Chiral phase transition for the flavor sector in this deconfined phase can happen and it has 
a beautiful geometric explanation: parallel profile of the flavor D8 and D8 branes stands for 
chiral restoring phase while connected U-shaped profile for chiral broken phase. In general, high 
temperature corresponds to chiral restoring phase while low temperature for chiral broken phase 
(the details on this topic can be found in fl7fl ). In the high temperature phase, the induced metric 
on the flavor branes has the following standard AdS form 

3 

dsj = (|) 1 (-f(u)dt 2 + d* 2 ) + 5 (j«L + . 2 ^ 2 ) . (5) 

The other model we are concerning here is the non-critical version of the above one. It is based 
on the supergravity description of N c coincident D4-branes in six dimensions with one dimension 
compactified on a circle as in the Sakai-Sugimoto model. The corresponding background geometry 
takes the following form in high temperature phase 

ds 2 = (|) 2 (-f(u)dt* + dz 2 + d X \) + (| ) 2 g- y (6) 

e = vm' r= t- /(!,) = 1 -(-)' t ~ (+ ' s = (+ ^ < 7 > 

Flavor quark sector can be introduced by adding Nf pairs of D4 and D4 branes into above back- 
ground geometry. One important feature of this model is that it does not have the undesired 
internal space, which may introduce the unwanted KK modes. Another striking characteristic of 
this model is that it has constant dilaton vacuum as for the D3-brane geometry in ten dimensions, 
which should signal that this model is a conformal one. As in the Sakai-Sugimoto model, chiral 
restoring phase means induced metric on the flavor branes takes the form 

R\ 2 du 2 



dsl = Q\-f(u)dt 2 + dx 2 ) + (^-J 



/(«) 



(8) 



For convenience in later calculations, we now rescale the radial coordinate u to make it dimen- 
sionless by a transformation ut/u = r. Then the induced geometry on the flavor branes takes the 
following simplified versions: 

dsl = (50 1 r-i (-f(r)dt 2 + dx 2 ) + flf ufr"! + r 2 <(S!j) , (9) 

e* = g s (jrfr-i, /(r) = l-r»; (10) 

*2 = © 2 i(-/W^ + ^) + (f) 2 ^, (ii) 

(12) 

Note that under the new coordinates, the interval for radial coordinate r is located in a finite 
regime: r G [0, 1] , and the horizon is now at r = 1 while boundary at r = 0, which will make later 
analysis convenient. 

Before concluding this subsection we give a short comment about one general feature of above 
models. Whether there is background on flavor branes, fluctuations of the flavor gauge field and 
scalar mode in chiral restoring phase will always decouple, which together with exact AdS forms 
of the induced metrics will greatly simplify later calculations. This is also one reason why we 
choose the transversely intersecting D-brane systems for study. As mentioned in Sec. |, the plasma 
structure function is completely encoded in the dynamics of the flavor U(l) gauge field propagating 
through the above-mentioned geometry, which is described by the DBI actions on the flavor branes. 

B. DIS: field theoretical definitions 

Deep inelastic scattering in QCD is a powerful tool in exploring the hadron structure. Here we 
mainly focus on the electromagnetic mediation between the charged lepton and the hadron. The 
basic objective of DIS is to compute the retarded current-current correlator defined by 

IV(fc) = i J d A xe- ik x e(xo)([Mx), J„(0)]>, (13) 

where k is the four-momentum of the electromagnetic current, (...) means quantum vacuum expec- 
tation, and Jfj,(x) is the mediated electromagnetic current. When considering the lepton scattering 
off the plasma, the hadron should be replaced by the plasma system and the vacuum polarization 
tensor Eq. (13) is modified to 

IV ( k > T ) = i J d 4 xe- ifc '^(x )([^(x), J„(0)]> T , (14) 



where {...)t means thermal expectation value in the plasma system. 
Though the gravity dual of the SU(N) x U(l) 

e.m. gsuge theory has not been established, we 
can use a non-dynamical electromagnetic field to model the photon just as in condensed matter 
physics given that electromagnetic coupling constant is very small. In the present context, the 
electromagnetic current is replaced by the flavor U(l) current and also denoted as Ju(x). Now 
we list the general structure of the thermal polarization tensor defined in Eq. (|b|). According to 
gauge symmetry and rotation symmetry of thermal field theory, it can be decomposed into two 
scalar functions as 

n M „ (k, t) = - ^\ n x (x, q 2 ) + u„ - kj^f\ k - kj^p) n 2 (x, Q 2 ) , (15) 

where rjuv = (—1,1,1,1), Q 2 = k^k^ is the current virtuality, and is the four- velocity of the 
plasma which will be chosen as = (—1,0,0,0) to signal that the plasma is at rest. We have 
also defined the Bjoken variable as x = —Q 2 /[2 (n ■ k)T\. Then the DIS structure function of the 
plasma can be extracted from the polarization tensor as 

F 1 = ^-lmIL 1 , F 2 ^ ~^ T k) lmU 2 . (16) 

For sake of later convenience, we now proceed to express the above two functions in terms of 
longitudinal and transverse polarization tensors TIll and H yy being introduced in Sec. Ill: 

F 1 = ^Imn ra = -*-Imn zz , (17) 

Z7T Z7T 

F 2 = ^ (^lmU LL + 2xF^j . (18) 



In obtaining these two expressions we have used the flavor current momentum defined in Eq. ( pop 
and assumed that the plasma is at rest. From Eq. ( |l8| ) we can see that in the interesting kinematic 
regime, to 2 /q 2 = 1 — Q 2 /q 2 ~ 1, F 2 can be simplified further to 

F 2 ~Q^l m n LL + 2xF 1 . (19) 

7T 



If the first term in Eq. (19) could be negligible, we then straightforwardly come to the familiar 
Callan-Gross relation F 2 ~ 2xF\. 

In particle physics, the structure function has been studied by the operator product expansion 
technique for specific hard processes. The parton model is suitable for weak coupling regime of 
high energy QCD, while in the present case for holographic quark-gluon plasma, which has been 
thought of as strongly-coupled one, we resort to its gravity dual for calculations of above quantities. 



There is a standard prescription [13| for calculating the retarded Green function such as Eq. (|l4]) 
under the approach of gauge/gravity duality. The current J M couples to its source (x, r = 0) as 

Sm = f d 4 xJ^(x,r = 0), (20) 

where An (x,r) is the flavor U(l) gauge field introduced in next subsection. The main idea under 
this prescription is to invert the operator Green function on field theory side to dual field Green 
function on gravity side, whose calculations just need classical gravity action. More details on this 
prescription can be found in [13]. Now we explicitly write down the expression for the polarization 



tensor defined in Eq. (14) in terms of the variables on dual gravity side 



a ^^ T ^ SA /Z M k) (2i) 



where S is the on-shell action defined in Eqs. (|26|) and (j34|) . Later, it will be found that due to 
coupling of A x and A t , we have to express the on-shell action in terms of longitudinal mode Al 



and transverse modes A y ,A z defined in next subsection. 



C. Basic equations: flavor current propagation in plasma 

As in many works on applications of gauge/gravity duality to strong coupling problems, we will 
choose radial gauge for gauge potential, i.e., A r = 0. This is enough for non-critical case while for 
the Sakai-Sugimoto model we also have internal symmetry on space. For brevity, we also set 
gauge potential along it to zero and assume that gauge potential does not depend on the internal 
coordinate. As to the action, we just retain to quadratic order in the gauge field fluctuation, which 
is enough for the propagation of the flavor current. In the following, we will just write down the 
main equations for the gauge fields (for details on the extraction of them, see [ |l4| , |l£j). 

For the Sakai-Sugimoto model, the DBI action for the fluctuation of the flavor U(l) gauge field 
takes the following form after integrating out the O4 space: 

S 8 = _l^£]_r 8iV/ y 54 J d A xdr v ^f f g MN g PQ F M pF NQ , (22) 

V~9eff = e _ V-dets50| 4 = g^R^r'* (23) 

where Fmp = OmAp — 8pAm is the field strength of the gauge field Am, g MN is the inversion 
of the induced metric in Eq. (||), the index M, N and so on just need to run the former five ones 
in Eq. (0), T$ is the D8 brane tension and V54 is the volume of unit sphere f^. The equation of 
motion (EOM) for the gauge field Am followed from this action is of Maxwell-type, 

9m (^g~JJg MN g PQ F NQ ) = 0. (24) 
9 



Now we turn to the momentum space by doing the following partial Fourier transformation of 
the gauge field 

d 4 k 



(x,r) 



e lkx A^k^r) , kf, = (-w, g,0,0) 



(25) 



(2vr) 4 

Without loss of generality, we have chosen the spatial momentum along just one spatial direction 
as done in the literature. Then in the partial momentum space, the on-shell action turns into the 
following form, 



27ri 2 ) 2 r d^k 

T 8 N f V si I ■ 7 ^ i y/^Jg rr {g tt A t (-k,r)d r A t {^ 



2 -"-■'•<>- J ( 2vr ) 4 
+ g*Ai (-k, r) 8 r Ai (k, r) } \ r r Z\ (i = x, y, z) , 



(26) 



while Eq. can be explicitly casted into following three ones 



uA' t + qf{r)A' x = 0, 

A ,, + dr(V^g77j9 rr 9 yy ) Xi _ g^_ i 



X '~!hJ]<r<i ini " g rr 



OJ' 



Ay = 0, 



2 9r {^glf f g rr g tt ) d r (^KJJg tt g xx 



yj-geffg rr g tt 



V-9effg u g 



tt n xx 



A'l + 



9? {V=g7J]g rr g u ) 



V=glfjg rr g tt 



a xx l 2 

g" 



" 2 \ d r ( v ^JJg n g--) d r {^Jfg rr g u ) \ A > =0 



V-g e ffg tt g xx ^-g e ffg rr g u 



(27) 
(28) 



(29) 



The component A z satisfies the same equation as the ^-component in Eq. (|~J), and we refer to 
them as transverse modes. Another important relation between A x and At is 



d r {V-9eff9 rr g u K) - V-9effg a g xx q M* + gA t ) = o. 



(30) 



In above equations, prime denotes derivative with respect to the radial coordinate r. Once the 
induced metric was inserted into Eqs. (28,29,30), they can be simplified to the following forms 

'Co 2 ~ q 2 f(r) f'(r) 



a" + 



A» + 



1 f'(r) 
2r f(r) 



a' + 



1 fir) 



2r f(r 
A L = ^q~ 1 rlf(r) (r^a)', 



r/ 2 (r) 2r/(r) 
_ q 2 f(r)-QJ 2 _ 
r p(r) v ~ 



a = 0, 



(31) 
(32) 
(33) 



where (w, q) 



4?rT 



(uj,q) are dimensionless variables, a denotes A' t , and Al is the longitudinal 



mode defined as A^ = qA t + u)A x . Note that we have written Eqs. (21 . . |,30) in a general form 
so that they should also be suitable for the non-critical model being studied later. 
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Similarly, we can write down similar equations for the non-critical version of the Sakai-Sugimoto 



model, but for brevity we just list the main results like Eqs. ( . . , ): 



4' + 



2 e<t> m J (2vr) 
- A t (-k, r) d r A t (k, r)] (i = x, y, z) 

f'{r) 1 



f(r) r 
f'(r) 1 



r(r) + &- em 



r 2 r J( r ) J2( r ) 



a = 0, 



_f(r) r 
4 1 „, , / 
R^q 



, Co 2 - q 2 f(r) _ 
A y H J2^ A v = °' 



(34) 
(35) 

(36) 

(37) 



Here the only difference from the Sakai-Sugimoto model in notations lies in the specific definitions 
of the dimensionless variables Co and q: 



CO 



CO 



(38) 



0.87rr' " 0.8vrT 

We now express the flavor brane actions in terms of Al,A v , A z and a as mentioned in last subsec- 
tion. Because similar calculations have been done in the literature for many times, we here merely 
list the final results. For the Sakai-Sugimoto model we have 

,2\2 



[ -^-T s N f V s4 J J^kV^Ff9 rr {^q- 1 A L {-k,r)a(k,r) 
+g u Ai (-k, r) d r Ai (k, r)} \ r r zl (i =y,z), 



(39) 



and for the non-critical model 



(2irl 2 s ) 2 T 4 N f 



T 



d A k 



2 R 3 J (2tt) a 

+f(r)A i (-k,r)d r A i (k,r)} \ r r zh 



r 1 {-q 1 A L {-k,r)a(k,r) 
(i = y,z). 



(40) 



Some remarks are due about these equations which will determine later calculations. Explicitly, 
these equations are more involved than the related ones in Ref. , which has already been confirmed 
in ||] for D3/D7 brane setup. We can attribute these complicatedness to the introduction of the 
flavor branes and the non-conformality of the model. However, once these EOMs were turned into 
standard Schrodinger types and taking high momentum and high frequency limits, we will find 
that these complications will automatically disappear and they are qualitatively in common with 
those of @ . 

We now proceed by following ]7j to turn these EOMs into standard Schrodinger types and 
discuss some general features of effective potentials respectively. We will find that these general 
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discussions for potentials will reveal some physical pictures for the DIS processes considered in this 
note. 

The equations obeyed by new fields have the following common form irrespective of the models 



(for time component a), 
: (for transverse component A y or A z ). 



(41) 
(42) 



For the Sakai-Sugimoto model, the explicit field transformations and effective potentials are as 
follows 



1 



Vl{t) 



A, 



V T {r) 



(1 — r 3 ) t 2 
1 



r (1 — r 3 ) 



3^ 2 



16r 



-3 - 78r 3 + 45r 6 ) + K 2 - q 



2 ^2 r 3 



T - 



1 — r 3 
1 



ijj T , 



r (1 — r 3 ) 



3^ 2 



16r 



(5 - 46r 3 + 5r 6 ) + if 2 - q 



2 r 3 



(43) 
(44) 

(45) 
(46) 



While for the non-critical model, we have following respective ones 



-4*1 



Vl{t) 

Ay = 

V T (r) 



(l-r 5 Y 



7-2 (-1 - 48r 5 + 24r 10 ) + K 2 - q 2 r 5 



(1 



1 _ r 



5^2 



4^2 



3 + 8r 10 - 36r 5 ) + K — q 



^ r 5 



(47) 
(48) 
(49) 
(50) 



In these equations, we have defined the dimensionless current virtuality as K 2 = q 2 — Co 2 . We will 
take this virtuality to be space-like, which amounts to saying that the process considered here is 
like the lepton deep inelastic scattering off the proton. But as noted in J?]], the final results will 
also be suitable for time-like virtuality if we take high momentum limit (especially if q 2 3> K 2 ). 
Since we are interested in the internal structure of holographic quark-gluon plasma, we should use 
high energy probe to explore this just as the DIS processes in pQCD. One basic difference between 
the plasma and a single hadron is that the former has an intrinsic scale (temperature). In short, 
we are focusing on the following kinematic parameter space: 



w»l, K 2 > 1 



(51) 
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Under the above kinematics, these effective potentials can be further approximated as follows- 7 ]. 
For the Sakai-Sugimoto model we have 



V L {r) -- 
V T (r) ~- 

while for the non-critical one 

V L (r) 

Vr(r) 



r (1 


— r 6 ) 






r (1 


_ r 3)2 




h- 1 


(1- 


- r 5 ) 2 




1 


(1- 


_ r 5)2 



16r 



+ K 2 - q 2 r 3 



16r 



4r 2 



7 2 r 5 



4^2 



+ if 2 - <f r 5 



(52) 
(53) 

(54) 
(55) 



Eqs. (0j42|,|52 

|)3],[54],^) are the main ingredients for later extraction of the structure functions of 
holographic plasma. 

We have seen that, in the interesting kinematic regime symbolled by Eq. (|5l|), for both the 
Sakai-Sugimoto model and its non-critical version, the effective potentials^] for longitudinal as well 
as transverse modes are qualitatively similar to those of the TV = 4 SYM case with(out) flavors 
0. [|. More explicitly, the maximum of longitudinal potential can be positive (corresponding to 
potential barrier), negative (with no barrier at all) or zero according to the value of q/K 4 (Sakai- 
Sugimoto model) or q 2 /K 7 (non-critical model). While effective potentials for the transverse modes 
are more involved because they start from a positive infinity and then falling to negative infinity 
very rapidly, which may complicate later analysis using WKB approximation to construct the 
wave function ipT- Recalling the fact that the dilaton vacuum for the Sakai-Sugimoto model is not 
a constant while for the non-critical version it is constant, these facts together seem to indicate 
that non-conformality of the Sakai-Sugimoto model is not essential on physical picture of high 
energy scattering process. This is a byproduct of the general behavior analysis for the effective 
potentials, later we will confirm this observation by direct extraction of the structure functions for 
both holographic models. 

Before concluding this section, we briefly summarize the physical picture in |?J governing high 
energy DIS from the viewpoint of non-relativistic quantum mechanics. For longitudinal mode, when 
the potential barrier builds in (corresponding to small spatial momentum case), the wave function 
cannot be imposed incoming wave condition at the horizon due to the high and narrow barrier, 

7 For the following four potentials, we have ignored polynomials of r in their final expressions. Strictly speaking, 
this is not right; but since we are in particular concerning their behavior near the boundary r — 0, which directly 
determines the polarization tensor, and their behavior at the horizon r — 1 is just for incoming wave boundary 

condition, we believe this is reasonable. 

8 We should have plotted these effective potentials for illustration as in FL H], but we skip it here for brevity since 

they are basically the same. 

13 



which indicating zero structure function. When taking into account non-perturbative tunneling 
effect, a small exponentially suppressed structure function can be obtained. We in following sections 
therefore merely focus on high spatial momentum limit. In this regime, the wave function will be 
complex one and incoming wave condition can be imposed at the horizon. Moreover, in this high 
energy regime, a partonic picture for the plasma exists. 



III. STRUCTURE FUNCTION OF HOLOGRAPHIC QUARK-GLUON PLASMA 



Now we have all the elements to calculate the polarization tensor defined in Eq. (|14|). As 
mentioned above, we should focus on the high momentum kinematics besides the one defined in 
Eq. ([y]). This means that the _£T 2 -terms in effective potentials can be ignored as well, which makes 
semi-analytical solutions for these Schrodinger equations possible. We in this section follow the 
standard WKB approach in non-relativistic quantum mechanics to construct these solutions. We 
present the calculations in detail for the Sakai-Sugimoto model and then list the final results for 
its non-critical version. 

We first discuss the longitudinal mode. Under high momentum approximation discussed in 
last paragraph, Schrodinger equation for longitudinal mode then takes following simple form near 
r ~ 0: 

^i(r) + + q 2 A ^(r) = 0. (56) 

Its general solution is a linear combination of the Bessel and Neumann functions: 

* L (r * 0) = j i/g (§pj + C2 ,~VV/ 2 iV 1/8 (^) , (57) 

and the constants c\, C2 will be determined by imposing incoming wave boundary condition at the 
horizon, which needs do matches of solutions at different regimes. 
Near r ~ 1, Schrodinger equation can be approximated as 



M(r) + . 2 Mr) = 0. (58) 

9 (1 — r) 



Its general solution is 



ipL (r ~ 1) = c 3 (1 - r)H§ + c 4 (1 - r)5 +i § . (59) 



Recalling Eqs. (43,45) and imposing incoming wave boundary condition at the horizon make us to 



conclude that C4 = 0, thus leaving the general solution near the horizon to be 

^ £ (r~l) = c 3 (l-r)5 _i i . (60) 
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Besides, we have to study in detail the solution suitable for intermediate regime far from the 
singularities at r = and r = 1. In this regime, Schrodinger equation is approximated as 

<M + TT^iMr) = 0. (61) 
(1 — r A ) 

For convenience, we now define the so-called canonical momentum p(r) and action s(r): 

qr f r f r qr 

P( r ) = q> s ( r ) = / p{r)dr = / dr ^. 

1 - r 3 J J 1 - H 

Then the two linear independent solutions in intermediate regime under WKB approximation are 

V4 1) (r) = ^=cos[s(r) + 1 ], (62) 

^\r) = ^=sm[s(r) + fo]. (63) 

The asymptotic behaviors for p{r) and s(r) at singularities are necessary for the solution matching 
underlying different regimes 

p(r ~ 0) ~ qr and p(r ~ 1) ~ - — -, (64) 

3(1 — r) 

~ 2 

or q 

s (r ~ 0) ~ and s (r ~ 1) ~ log (1 — r) + constant. (65) 

2 3 

The next step is to match these solutions to determine c\, c 2 and C3. In doing this, we need the 
asymptotic expansion for the Bessel or Neumann function for very large variable^ 

Jl " (ir) = cos (x " Te ~ l) ' (66) 
(f -5-5) ■ «-» 

It is then easily found that if we choose 4>i = <p2 = — 5vr/16 and c 2 = ici, the matching of solutions in 
different regimes accomplishes. Note that the condition c 2 = ic± is the main result from the solution 
matching, which is also a direct reflection of the incoming wave boundary condition imposed at 
the horizon. 

Now we easily arrive at the boundary behavior for the wave function V , i( r ) as 
i, L (r,0) = Cl ^/4 r i/2 Ji/g ^ + ^^1/2^ ^ 
,1/4,1/2^(1) fqr^_ 



cl?WH±% (68) 



9 Although r is small here, q is large and we therefore regard qr 2 /2 to be large enough. 
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In the second line of above equation, we have written the solution as the first kind Hankel function 
with order 1/8. 

The effective potential for transverse mode is more involved, but the analysis under WKB 
approximation is similar, therefore we here just list the final results. Matching of solutions in three 
different regimes (near horizon, near boundary and intermediate regime far from these singularities) 
results in the following boundary behavior for transverse mode A y and A z 



V> T (r * 0) = cl^y (gr 2 /2) • (69) 

The only undetermined constants c\ and c\ can be expressed in terms of boundary values of 
the gauge field Al{t = 0) = Al(0) or A V)Z (r = 0) = Ai(0) respectively. This can be achieved by 



using Eqs. ( . J |) and the final results essential for later extraction of structure functions are 



listed as below. For longitudinal mode we have 



Mr - 0) = Cl grtW/*(l - ,3) J^^H® (f ) \, (70) 



«(-0)-iW7f^<f?l (71) 



(1-r 3 ) V8 V 2 

Cl - Wm^^ ALm (72) 

and for transverse modes Aj (i = y,z) we have 



Mr^o)=^^L H m(sfj, (73 , 

* = ¥^m M0 >- (74) 

With these solutions, we can derive the expressions for the thermal polarization tensor defined in 
Eq. flnj): 

^(k,T) = T ^^ XN f N c T^, (75) 

lmU yy (k, T) = Imn 2Z (k, T) = —^L_XN f N c T 3 f/\ (76) 

and other components are exactly vanishing. Referring to above results for the polarization tensor, 
we have just listed the imaginary parts which are directly related to the structure function. In fact, 
the real part of polarization tensor are divergent and therefore needs regularized (here we skipped 
these details). Moreover, we here for simplification assume the regularization will not introduce 
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new terms into the imaginary parts. Therefore, the structure functions F\ and F 2 are easily derived 
via Eqs. dT^H,^) 



F 2 (k,T) 



uj 2 



,2 



1/4 

7T i2r 2 (i/8) 5s ' 



2xF l . (78) 



In the first line of Eq. (f7l) one can easily show that the first term can be ignored when comparing 
to the second one in the interesting kinematic regime q/K* > 1. Following §, we now express 
the two structure functions in terms of the Bjoken variable x defined in Sec. | 3 and the flavor 
current virtuality Q 2 as 

Fl ( Sj q2) = \NfN c T 3 ( -^^] ~ XN f N C T 3 ( (79) 

v ^ ' 108r 2 (3/8)g s ; \8ttxT 2 J 1 \8itxT 2 J ' v ' 

f SO 2 \ 3/4 

F 2 (x,Q 2 )~2xF 1 (x,Q 2 )~2\N f N c T 3 x( ■ (80) 

In obtaining the above two equations, we have used one approximate relation Co ~ q to express the 
spatial momentum as q ~ Q 2 /(2xT). 

Before closing this subsection, we briefly carry out similar analysis for the non-critical model. 
Since in the previous subsection the general procedure for using WKB method to solve Schrodinger 
problem has been presented in detail, we here just list corresponding key results. The solution for 
longitudinal mode near boundary r = behaves as 

Mr - 0) = d^-y/Ml - r 5 ) U^H® (^> 7/2 ) } , (81) 
a(r ~ 0) = Cl q^-^==H^ (^qr^ , (82) 
<?i = -y^^4x(0); (83) 



and for transverse mode Ai(r) (i = y, z) 



(84) 



ci = r£§jT) A ^ (85) 

The polarization tensor can be derived by inserting these solutions into the on-shell action 
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Eq. (EC 



lmU LL {k,T) = 1 (86) 

^ yy (k,T) = 25T2{2/7W T . (87) 



Then the structure functions are 



Q 2 x V3N f N c | jVERNfNq ^ ( q\ 4/7 ' 



uj 2 

~2xF 1 (k,T). (89) 

In the first line of Eq. (Q]), one can easily show that the first term can be neglected in the 
interesting regime (? 4 / 7 3> K 2 , so we have the approximate equality of the second line as in the 
Sakai-Sugimoto model. We now express these results in terms of the Bjoken variable x and the 
flavor current virtuality Q 2 as in the Sakai-Sugimoto model 



( ^o 2 \ 4/7 

F 2 (x,Q 2 )^2xF 1 (x,Q 2 )~2N f N c T 2 xl ^-^J • (91) 



Eqs. (["3 "3 ^Ol"HOO) comprise the main results of this note. Now we have a short 

remark on these results and also a brief comparison between these two models as well as the 
well-investigated N = 4 SYM case. 

The first point is that we have an analogy of Callan-Gross relation F 2 — 2xF\ in the interesting 
kinematic regime considered in this note. In pQCD, this relation holds only at relatively large 
Bjoken variable x, where parton structures of hadrons are dominated by the valence quarks. This 
relation has been already obtained in DIS off W = 4 SYM plasma with(out) flavors and here 
it also holds for the Sakai-Sugimoto model as well as its non-critical version. Since the setups of 
holographic dual of sQGP are very general, and what is more is that the physical picture underlying 
the Tt or the flavor current DIS off the plasma is quite simple and general, these surprising facts 
seem to indicate that it may be a general relation for holographic quark-gluon plasma. 

The second key point is concerned with the non-conformal characteristic of the Sakai-Sugimoto 
model, which is also the main motivation of the present study. Because in our interesting kinematic 
regime, the two structure functions are related to each other by the Callan-Gross relation, we then 
mainly focus on the function F\. It is clear that for two models F\ presents scaling behavior on 
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their dependence on temperature T and dimensionless spatial momentum q. Their dependence on 
dimensionless momentum q are approximately the same, while on temperature are quite different: 
~ T 3 for the Sakai-Sugimoto model and ~ T 2 for the non-critical one. Recalling that the latter 
scaling behavior ~ T 2 has also been valid in |?J for the M = 4 SYM plasma, we may think of the 
essential effect of non-conformality of the Sakai-Sugimoto model as the T 3 -scaling behavior of the 
Fi structure function. However, this guess needs further confirmation or cancelation because the 
non-conformality of the Sakai-Sugimoto model is not well-controlled. Moreover, the gauge coupling 
constant of strong interaction has a logarithmical running with evolving energy scale and simple 
gravity realization of this kind of gauge theory has been established by combining top-down and 
bottom-up approaches in it may be interesting to resort to this kind of model to probe the 
effect of gauge coupling running on the internal structure of sQGP. 

The last point we would like to stress is about the pre-factors for the structure functions. 
Similar to S, NfN c counts the number of freedom of the plasma, and we here probe the quark 
sector. The models considered in this note has left some stringy imprints on the quantities of 
field theoretical side (here they are the plasma structure function). This can be easily read from 
Eqs. Q)- which have explicit dependence on string coupling constant g s and string length 

l s . Moreover, the behavior concerned to these two stringy parameters seems to be different between 
the Sakai-Sugimoto model and its non-critical version. These facts seem to be very strange and 
even unaccepted because we here focused on the field side quantities and they should not show 
explicit dependence on the parameters on the gravity or string side. Compared to related results 
of the W = 4 SYM plasma, these undesirable behaviors have not come out there, which seems 
to say that the D3-brane geometry is a much better gravity dual of some QFT in describing field 
theoretical physical quantities under gauge/gravity duality technique. However, if we recall that 
the models we take here are different from realistic QCD theory, then it is acceptable that the 
results are counter-intuitive from the viewpoint of field theoretical considerations. 



IV. SUMMARY AND OUTLOOK 

In this note, we have used high temperature version of the Sakai-Sugimoto model, a quite 
successful gravity dual model of QCD-like theory, to explore the internal structure of holographic 
quark-gluon plasma. The physical process we analyze here is like the well-investigated DIS in 
standard QCD, but with the scattered proton replaced by the plasma system and the mediated 
electromagnetic current simulated by the flavor current as in Q . We have seen that the procedure 
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under gauge/gravity duality technique to study DIS off holographic quark-gluon plasma is quite 
general and easily promoted to other gauge/gravity duality setups. The result obtained in this 
note for the structure function under the Sakai-Sugimoto model is quite different from the well- 
studied J\f = 4 SYM plasma with(out) flavors. This should be regarded as the effect of the 
non-conformality of the Sakai-Sugimoto model, which is the most important motivation for our 
present study. To confirm this, we have also chosen the non-critical version of the Sakai-Sugimoto 
model for a comparative study. We found that the structure functions for the latter model and the 
N = 4 SYM plasma are much alike. The result of this note seems to contradict the intuition from 
the fact of asymptotic freedom of pQCD. But we should keep in mind that holographic quark-gluon 
plasma considered in this note is a strongly-coupled system and we should not expect it to behave 
exactly as the weakly-interacting regime of realistic QCD. In fact, a more realistic holographic 
QCD model taking care of the running of the gauge coupling constant was proposed in | |19| and we 
expect to use this model to explore the effect of gauge coupling running on the structure of sQGP. 



In realistic QCD, Nf and N c are of 0(1), while the applicability of gravity dual of SU(N C ) 
gauge theory usually requires a large N c limit. Therefore, we cannot directly compare our results 
with the data from heavy-ion collision. One prescription to overcome this obstacle is to consider 
the flavor backreaction to the background geometry and then carry out similar calculations in 
the Nf/N c ~ 1 limit. Although the hadronic matter produced in heavy-ion collision is of high 
temperature and high density, we in this note only take into account the high temperature element 
as in the literature. So in this sense, our present models are not so realistic and need to be promoted 
to high energy and high density quark-gluon plasma case. Fortunately, under the gauge/gravity 
duality setup, the matter density also has a gravity realization — 7Z charge or flavor charge, which 
can be realized by rotating color brane along the internal space or as the time component of the 
flavor gauge field respectively. Then, the analysis of DIS off quark-gluon plasma at high density 
and high temperature can also be carried out by including this element in the model setup. We 
leave these problems for future investigations. 
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